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Patent Application 



A COMMUNICATION SYSTEM AND METHOD 
FOR PERFORMING FAST SYMBOL ESTIMATION 
FOR MULTIPLE ACCESS DISPERSE CHANNELS 



The present invention relates to wireless and wireline communication systems. More specif- 
ically, the present invention relates to an improved method for estimating symbol sequences 
of multiple sources sharing the same communication media. The invention is particularly 
useful for code-division-multiple-access (CD^LA.) receivers where the channel length is only 
limited to two neighboring symbols. 

In a scenario where m receivers are utilized, the received signal at each receiver is the simi 
of the output of channels driven by u different user signals plus the noise at the receiver. The 
channel between receiver i and user j can be denoted dshij. In general, hij # /i^ if i 7^ g or 
j ^ n. Noise signals at different receivers are assumed to be independent of one another. If 
the user signals are actually CDIVLV signals and the largest duration of /i^ , 2 = 1, • • ■ , m and 
J = 1, * • • , It is short enough compared to the symbol duration of those CDMA signals, the 
orthogonality of the spreading codes can guarantee the separation of CDMA signals at each 
receiver. Intersymbol interference (ISI) and multiple access interference (MAI) are negligible. 
When the largest duration of the signature waveforms is comparable with, or even larger 
than the symbol duration, ISI and MAI cem become severe. The performance of detection 
made on symbols individually will not be good due to ISI and MAI. Such a situation calls 
for joint detection by which channel effects are modeled adequately and symbol sequences 
are estimated jointly. In minimum-mean-square-error (MMSE) based joint detection, which 
will be referred to as joint detection in the following, an inverse filtering process is used to 
remove intersymbol interference and multiple access interference. To do this, a coefficient 
matrix H formed from the channel responses is inverted. More precisely, the pseudo-inverse 
of this coefficient matrix: {H^H)~^E^ must be found. The main problem associated with 
joint detection is that the computation thus involved is intensive. The dimension of if is 
determined by the number of users, the length of spreading sequences, the length of symbol 
sequences, and the number of receivers involved (even with only one receiver, the effect of 
multiple receivers can be achieved by sampling a received continuous signal at a multiple 
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times of the symbol rate.). With a moderate number of users and a moderate number of 
symbols, the dimension of the coefficient matrix can be large. To speed up the inversion 
of the coefficient matrix, a Cholesky decomposition on H*H qbh be used. However, the 
computation amount involved in the Cholesky decomposition can still be substantial, as the 
requirement on matrix operations such as matrix inverses increases linearly with the number 
of symbols. It is therefore desirable to have a method to speed up the joint detection 
computation. 

Prior to this invention, a class of multiuser detectors have been developed. The most 
prominent ones among many others include It Lupas and Verdu, "Linear multiuser detec- 
tors for synchronous CDMA channels", IEEE Trans, on Information Theory, 1(35): 123-136, 
January 1989.; Z. Xie et al. "A family of sub-optimum detectors for coherent multiuser 
communications", EEEE Journal of Selected Areas, in Communications, pages 683-690, May 
1990; and Z. Zvonar and D. Brady, "Suboptimum multiuser detector for syncrhonous CDMA 
frequency-selective Rayleigh fading channels", IEEE Trans, on Communications, Vol. 43 
No. 2/3/4, pp.154-157, February/March/ April, 1995, Angel M. Bravo, "Limited linear can- 
cellation of multiuser interference in DS/CD^LV asynchronous systems'', pp. 1435-1443, 
November, 1997. They either design various finite-impulse-response filters to process the re- 
ceived signals, which are in essence approximate solutions of the MMSE based joint detection 
solution, or use Cholesky decomposition to solve the equations involving the large coeflScient 
matrix, e.g. Paul D. Alexander and Lars Rasmussen, *^n the windowed Cholesky factoriza- 
tion of the time- varying asynchronous CDMA channeT, IEEE Trans, on communications, 
vol. 46, no. 6, pp.735-737, June, 1998, which gives the exact MMSE joint detection result 
with high computational complexity, i.e. complexity of matrix operations such as matrix 
inversion increases linearly with the nimiber of symbols. 

As in many communication systems, the length of ISI and MAI is often limited to sev- 
eral symbols, besides being Toeptliz or block Toeplitz matrix, the coefficient matrix is also 
banded. Fast algorithms developed for generic Toeplitz or block Toeplitz matrices, e.g. 
Georg Heinig and Karla Rost, "Algebraic Methods for ToepUtz-like Matrices and Opera- 
tors", Birkhauser, 1984, do not take advantage of the fact that the coefficient matrix is 
banded. The afore-mentioned Cholesky factorization based algorithm, while taking advan- 
tage of the banded nature of the coefficient matrix, does not take advantage of the fact 
that the coefficient matrix is also Toeplitz or block Toeplitz, which leads to considerable in- 
crease of computational demand as the number of matrix operations such as matrix inversion 
increases with the number of symbols to be estimated. Simpler versions of the method dis- 
closed in this invention, have been used to solve numerically Poissdn equations, e.g. Roland 
A. Sweet, "A cyclic reduction algorithm for solving block tridiagonal systems of arbitrary 
dimension", SIAM journal on Numerical Analysis, Vol. 14, pp. 706-720, September, 1977. 

SUMMARY OF THE INVENTION 

The present invention discloses a symbol detection system and method for multiple users 
in scenarios where the medium is shared by multiple users. The symbol detection system 
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and method involve fonnulating the received signals as the product of the inverse of a block 
triadiagonal matrix and a vector formed by the symbol sequences. The block tridiagonal 
matrix is derived from the channel responses. The present invention includes a fast algorithm 
to find the product. The computational complexity of the disclosed method is linear with 
the number of symbols and the requirement for matiix operations such as matrix inversion 
increases with the logarithm of the number of symbols. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a graphic illustration of a wireless communications setup for the disclosed technique. 

FIG. 2 is a graphic illustration of a wireline communications setup for the disclosed technique. 

FIG. 3 is an overview of the disclosed algorithm. 

FIG. 4 is a flow chart of the forward reduction procedure. 

FIG. 5 is a flow chart of the backward substitution procedure. 

While the invention is susceptible to various modifications and alternative forms, spe- 
cific embodiments thereof are shown by way of example in the drawings and will herein 
be described in detail. It should be understood, however, that the drawings and detailed 
description thereto are not intended to limit the invention to the particular form disclosed, 
but on the contrary, the intention is to cover all modifications, equivalents and alternatives 
falling within the spirit and scope of the present invention as defined by the appended claims. 

DETAILED DESCRIPTION OF THE DRAWINGS 
INTRODUCTION 

The basic problem under consideration is that of estimation of symbol sequences with inter- 
symbol interference and multiple access interference. To give a concrete setup of the joint 
detection problem under study, we consider a wireless CDMA system with multiple anten- 
nas (receivers) at the base station, and consider the jomt detection problem of finding the 
symbol sequences from these users, i.e., the joint detection problem of the uplink signals 
received at the base station. As was made clear previously, by oversampling, the effect of 
multiple virtual receivers can be achieved with only one physical receiver, the applicability 
of the procedures which will be given is apparent for the joint detection problem of downlink 
signals also. In a code division multiple access system with J users, there are M {M > 1) 
antennas and hence M physical receivers at the base station, and the length of the spreading 
sequences is P. The received signal at each physical receiver is oversampled at Mt times the 
chip rate thereby producing a plurality of received signals each sampled at the chip rate. 
In other words, for each physical receiver, the received signal is oversampled to produce a 
number of chip rate signals. Hence at the base station, the signal received at antenna i can 
be represented as follows: 
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where hik{n) is the signature waveform of the channel between the A;-th user and the i- 
th antenna at the base station. The signature waveform can be identified either by blind 
methods or with training sequences, e-g. Tong et al., "Blind Identification and Equalization 
Based on Second-Order Statistics: A Tmie-Domain Approach'', IEEE Trans, on Information 
Theory, March, 1994, 

FIG. 1 shows one example of an embodiment where the invention may be used, A base 
station (102) receives signals from wirdess terminals (104, 106, 108, 110). For some wire- 
less terminals(106, 108), there is more than one path between them and and the base station, 
as made clear by the reflective paths by reflectors (112, 114). Path length difference among 
multiple paths can destroy the orthop)nality of signal waveforms of different transmitters 
which may have existed. Methods based matched filters and simple despreading are not able 
to offer satisfactory results under such circumstances. MMSE joint detection can offer much 
better results with heavier computation complexity. The disclosed method is able to reduce 
the computation amount requirement for implementing MMSE joint detection considerably. 

In the embodiment of FIG. 1, the symbol estimator of the invention may be comprised in 
the base station (102) for estimating symbols in uplink transmissions from wireless terminals 
(104,106,108,110). The symbol estimator of the invention may also be comprised in one 
or more of the wireless terminals (104,106,108,110) for estimating symbols in downlink 
transmission from base station (102). 

FIG. 2 shows another embodiment of the disclosed invention. A cable modem node 
(200) receives signals from multiple cable modems (206, 208, 210) which share the same 
cable medium (212) to communicate with the cable modem node (200). The channels from 
different cable modems to the cable modem node (200) will be different. In the embodiment 
of FIG. 2, the symbol estimator of the invention may be comprised in the cable modem node 
(200). Similar to the situation shown in FIG, 1, the disclosed technique can also be applied, 
to speed up the joint detection method that combats the dispersive channel effects. It is not 
diflBlcult to see that the invented technique can be applied in a reversed way, i.e., the cable 
modems (206,208, 210) may include the symbol estimator of the invention and may receive 
the signal from the cable modem node (200). 

FIG. 3 shows one embodiment of the method used by this invention to estimate symbol 
sequences of multiple users. The base-band signals are collected and sampled from multiple 
receivers (302). In (304), At each receiver, the channel response between a transmitter 
and that receiver is used to filter the sampled received signal at that receiver. As there are 
J channel responses between all the txansmitters and that receiver, there are a total of J 
filtered outputs produced from the sampled received signal at receiver. All the outputs are 
multiplexed to form a single data vector, according the transmitter indices, in an ascending 
order. Next all the multiplexed data vectors from all the receivers are added up to form one 
single vector Yi (318). 

Channel responses hmj{n), m = 1, M, j = I, • • • , J (312) between transmitters 
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and receivers are identified either by training sequences or other methods which are not 
specified in this invention. Spac^time correlation matrices T;, z = ci, d — 1, - • * , — (d — 
1), — d, are constructed firom those channel responses (312). Channel description matrices 
Ai, Bi, and are block Toeplitz matrices, the blocks of which are defined by T^, i = 
d, d - 1, • • • , -{d - 1), -d, which are constructed to form the description matrices (316), 
Forward reduction steps (306) are initialized with the second set of matrices and vectors, i.e, 
Ai, Bi, Di and Yi. In the course of forward reduction, using the bisection approach, some 
intermediate matrices which will be used in the badcward substitution steps in (308) are 
stored. At the beginning of backward substitution steps (308), the intermediate solution 
is foxmd. Solved unknowns and intermediate matrices which have been stored in the forward 
reduction steps are used to find otha^ unknowns. The output of the backward substitution 
steps (308) is demultiplexed (310) to produce the estimate of symbol sequences sent by 
transmitters. 

FIG. 4 is a flow chart showing the forward reduction procedure. 
FIG. 5 is a flow chart showing the backward substitution procedure. 



DATA FORMATION 

In this disclosure, cases where signature waveforms of arbitrary duration are considered. In 
the following, we will write the received signals as the convolution of signature waveforms 
and oversampled symbol sequences in a structiiral way, which facilitates the derivation of 
the fast algorithm. If the signature waveform duration is limited to + 1 symbol duration, 
then we can pad zeros to hi^k{n) so that the length of /it^(n) is precisely {d + 1)P. Assume 
that hmj{n), 1 < m < M, 1 < j < / are column vectors, which are identified with either 
training sequences or other means, which is not specified in this invention (312). We first 
define some notation. We divide hmj{n) into {d -h 1) parts: 



Knj = 



where h^^^\ j = l, d + l, arePxl vectors. We put together the signature waveform 



vectors: 



(mJ) 



H(ni) 
d+1 
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where 



(ml) 



(mJ) 



h£°^ is a P X J matrix, k = l, d + l. 

Put the received signal at receiver m into a vector fashion: 



,r£(i + d)P) J 




xx(l) 
:«:2(1) 

xj{l) 

X2{L) 



By stacking the received signals: Rm, m = l, • • • , M, one over another, we have 



\Ri 1 








• Wi ■ 


R2 




H2 


x + 


W2 








✓ V 




' V 

R 




H 




w 



The objective of joint detection is to find X. If 

{L + 1)PM>LJ, 

then nonnally H is a full rank matrix. 
The zero forcing solution is given by 

Af M 
m=l m=l 

The MMSE (minimum mean square error) solution is given by 

M M 
m=l m=l 



where a- is the variance of the noise. We define a matrix: 



M 



and a vector 



M 



where v can be either 0 or cr^ depending on whether the zero-forcing solution or the MMSE 
solution is in need. B^Rm, m = 1, • • - , M are formed in (304). Y is formed in (318). In 
both zero-forcing and MMSE solutions, the computation ofX can be separated in two steps: 

1. Find Y. 

2. FindT-iy. 

In either case, T is a block banded matrix, defined in the following way, 

To 



T = 



(5) 



where 



° JLJxLJ 



Ti = E E(HSl-*)'Hi':L,^ + 5{{)vl, T.i = i = 0, -1 • • • , -d. 

m=lk=Q 



And Ti, i = d, d - 1, • • • , -(d - 1), —d are called space-time correlation matrices (314), 
delta{-) is the Kronecker function. 

We can represent T as a block tridiagonal matrix, if d divides L: 



1. the diagonal block A is given by 

To 

A= ■ 



T_M_ 



(d-i) 



dJxdJ 
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2. the supra-diagonal block B is given by 



B = 



■idJxdJ 



We can write 



TX = Y 



(6) 



as 



A B 
B* A B 



B* A B 
B* A 



X = Y. 



(7) 



FAST JOINT DETECTOR 
DERIVATION 

The system of equations as given in (7) can be treated as the special case of a more general 
formulation defined by a triplet [A, B, D]: 

Li 



B-i A, Bi 



B-i Ai Bi 
B'i Bi 



Xi = Yi, 



(8) 



where A,, B,, and are J x J matrices, which are called channel description matrices 
(316), 



Xi = 



(9) 
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Yi,i 



(10) 



with Xij and Yij being J x 1 vectors. 

We linearly combine equations for the following cases: 



1. for 2k = 2 (i.e. Jk = 1), 

2k -1: -B'iA"^ 
2k: I 



X 
X 



/r 



2A; + 1 : -BiA; 



-1 



A, Bi 
B'i Ai Bi 

B'i Ai Bi 



Xifik-l 
Xi,2k 
Xi,2k+l 



Yifik-l 
Yi;2k 



(11) 



and we have 

[O -B%A-^Bi + Ai-BiA-^B*i 0 -BiAT^B^ ] 

= -B%vy,^-i + yi,2Jk-B,Ar'yi,2/:+i, 



(12) 



or 



(13) 



[ -B*iA-^Bi + Ai- BiA-^B-i -B^ A'^B^ ] 
= -B'iA-^Yi^-i + Yi,2k-BiA-'Yi,2k+i. 

2. For every even index 2A;(4 < 2ifc < Ij - 3), we can linearly combine the 2k - 1**, 2A;*'*, 
and 2k + 1** equations in system of equations (8): 



2A: - 1 : -B'iA"^ 
2k: I 
2A: + 1 : -BiA"^ 



and we have 



X 
X 



B'i Ai Bi 

B'i Ai Bi 

B'i Ai Bi 







\ 


















Xi^zk+i 















(14) 



[-B^A^'B^ O -BvA-^Bi + Ai-BiA-^B% 0 -B^A-^Bi ] 

= -B%Ar'yi^.i+yi,2jfe-B,Ar'n2Jt+i, 



^t,2Jfe-2 
^i.2jk-l 
^i,2Jt 
-^i,2A:+l 
_ -^i,2Jfc+2 



(15) 
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or 



(16) 



3. Depending on whether Li is odd or even, we need to study two different cases. 



(a) If ii is odd, the last even index 2k between 1 and Li is Li- 1. Let Li = 2j + 1, 
so ik = j. We combine the last three equations 



-1 



2j-l: -B*iA-^ X 
2j: I X 

2j + l: -BiDi 



we have 



/r 



B'i Ai Bi 

B'i Ai Bi 
B*i Bi 



Xi^j-2 

Xi^j+1 



(17) 



or 



[-B%Ar'B% 0 -B^A^'Bi + Ai-BiD^'B^ o] 

[ -B*iA-^B% -B*iA-^B^ + Ai - BiD-^B'i ] 
= -B*iA^^Yi^j^i + Yi^j - B*iD~^Fi,2j+i- 



(18) 



(19) 



(b) If Li is even, the last even index 2k between 1 and is Li. Let Li = 2j. We 
combine the last two equations: 



2j-l: -B\A-^ X 
2j: I X 



B'i Ai Bi 
B'i Di 



Xi,2j-l 



, (20) 



we have 

[-B-iVB'i 0 -B-iA-^Bi + Di ] 



■X"»,2i-2 



= -B'iA-'Yi,2j-i + Yi,2i, 

(21) 
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or 



= -B^A^'n2i-l+l^^^;. (22) 



To facilitate a structural description of the aewly-generated equations, we can define the 
following: 

Ai+i := -B%A-^B, + A,-BiAr'B%, 
B^+i ~ -B,A-'Bi. 

If Li is odd, we define 

Di+i - A, - BiDr'B% - B%A-^Bi, 

and 



If Li is even, we define 
and 

We also define 



I,+i:=-(ii-l). 



(23) 

(24) 
(25) 
(26) 
(27) 
(28) 



Xi^i^ := Xi^j 1 < 2k < Li, 

At+i, Bi+i, and Dj+i are the new channel description matrices. 

Substitute the above newly-defined variables into equations (16) and (13) for 1 < 2A; < 
Li - 2, and equation (19) if Li is odd, or equation (22) if Li is even. We obtain the following 
system of equations: 



At+1 Bi+1 
B*i+i Ai+i 



B 



B*i+i Ai+i Bi+i 
B*i+i Di+i 



= Yi 



(29) 



We call the system of equations given in (29) as the i + system of equations. Hence a 
system of equations of reduced dimension involving only Xi^ of even indices as unknowns is 
obtained in the cancellation process, and Yi^i is obtained. It can be seen that the coefficient 
matrix defined by the triplet [Ai+i, Bi^u Di+i] is still a block-tridiagonal matrix, even 
though of smaller dimension than the coefficient matrix defined by by the triplet [At, B^, Di]. 
Hence the same procedure can be again applied to this new system of equations. With the 
execution of the procedure each time, the size of the system of equations is reduced by at 
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least a half. If we start from the system of equations as given in (10), by executing 
this procedure not more than logziLi) times, we can generate the 2^^, 3'''', - • • systems of 
equations. Eventually, at some 5, the s-l^ system of equations would be given either as 



B*,.i A^_i Ba-i 
BVi D.-iJ 



or 



A,_i Ba-i 



(30) 



(31) 



(30) and (31) can be reduced to the ^ system of equations involving one u x 1 subvector: 
Xs,i as unknown, respectively tiurough (17) and (20): 



BsXs = 



we have 



(32) 



(33) 



The just-found X,,i(or X^-i^ as in the s - 1** system of equations) can be substituted into 
the s - 1^ system of equations. In general, for the i** system of equations, from solution 
of the i + 1^* system of equations, Xij of even indices are already known, and Xij of odd 
indices are still unknown. We have three cases to consider. 



1. Consider the first equation in the t*'' system of equations: 



(34) 



where Xi^{Xi+i,i) has been found from the solution of the i + 1'* system of equations. 
We have 

Xa = ArHyi.i-Bi,i^m.i)- (35) 



2. Consider the 2k + 1^* equation in the system of equations, A; = 1, 2, etc. 



2A: + 1 : [ B*,- A, B^ ] 



(36) 



where the values of Xi;2t{Xi+iji) and Xi,2k-i-2{Xi+i,k+i) have been found from the solu- 
tion of the i + V* system of equations. We have 



or 



(37) 
(38) 
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3. If Lj is odd(Li = 2k + l), then Xij,, is found from the last block equation: 



[B-i Di] 




(39) 



with 



(40) 



Consequently for the system of equations, we just solve the block equations corre- 
sponding to odd in<fices. Hence all the unlmowns in the system of equations are found. 
The above procedure can be recursively appBed to the i — 1^, i — 2"^, • • • sjrstems of equa- 
tions until the unknowns in the system of equations as given in (7) axe all found. 

Apparently, the main computation is on tiie right-side of the equations. The update of 
the triplets involves li^t computation. 

In summary, the solution of the block tridiagonal system of equations as given in (10) is 
executed with two procedures. In the first procedure, the size of the system of equations is 
reduced, which is called forward reduction (306), a flow chart of which is shown in FIG. 4. 
In the second procedure, solved unknowns are substituted into larger systems of equations, 
and more unknowns are solved, which is called backward substitution (308), a flow chart of 
which is shown in FIG. 5. 



The whole procedure is naturally separated into two parts: the forward reduction procedure 
and the backward substitution procedure. 
The forward reduction procedure: 



1. Ai := A, Bi ~ B, Di := A, Yi := Y, Li := i := 1. 

2. Find the inverse of A^: A,"^. 

3. Qi := BiA-^ and := B^iAT^ 

4. Ai^i :=A,-EiB,-(E,B,r. 

5. Bi+i := — QiBi. 

6. • If Li is odd, 



COMPUTATIONAL PROCEDURE 



(a) P,:=B,Dr\ 

(b) D,+i-A,-PiB%-E,Bi, 

(c) 1,^1 := 

(d) := —EiYi^^i + Yi^2k — Qi^;2Jfc+ij 1 < ^ < Li^i, 

(e) Yi^ij^, := -E,yu.-2 + Yi,L,., - FiYi,^,,. 



• If Li is even, 
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(a) Di+i :=Di-E,Bi, 

(b) Li+i := ^, 

(c) Yi+i^ := -EiYi^.i + ^2* - Qili^+i, 1 < A; < Xi+i, 

(d) := -Ei5^irt<-i + Yi^- 

7. Store A,:^\, D,:^V Qi. Et> P» (if is calculated), and Yi. 

8. i := i + 1. 

9. Repeat 2 to 8 until Li = 1. 
10. s:=z. 

A flow chart of the above-stated procedure is given in FIG. 4. 
The backward substitution procedure: 

1. Solve DsXs = Y,. 

2. i:=3-l. 

3. Xi^2k := Xi+i^, 1 <k < Li+i. 

4. j := 1. 



5. := < 



f A-%-ErXij+i, ifj = l. 

-Q:Xij_i + A-^y^j - E^Xij+i, if K J < ii. 
I -P?Xij_i+D-^yij, if j=Li. 



6. j := J +2, if j < ii go to 5. 

7. i := i - 1, if i > 0 go to 3. 

8. is then the estimate of X. 

A flow chart of the above-stated procedure is given in FIG. 5. 

Although the system and method of the present invention has been described in connec- 
tion with the preferred embodiment, it is not intended to be limited to the specific form set 
forth herein, but on the contrary, it is intended to cover such alternatives, modifications, 
and equivalents, as can be reasonably included within the spirit and scope of the invention 
as defined by the appended claims. 
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